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Abstract 

In this paper we compute the scaling functions of the effective central charges for vari- 
ous quantum integrable models in a deep ultraviolet region R — > using two independent 
methods. One is based on the "reflection amplitudes" of the (super-) Liouville field theory 
where the scaling functions are given by the conjugate momentum to the zero-modes. 
The conjugate momentum is quantized for the sinh-Gordon, the Bullough-Dodd, and the 
super sinh-Gordon models where the quantization conditions depend on the size R of the 
system and the reflection amplitudes. The other method is to solve the standard ther- 
modynamic Bethe ansatz (TBA) equations for the integrable models in a perturbative 
series of 1/ (const. — lni?). The constant factor which is not fixed in the lowest order 
computations can be identified only when we compare the higher order corrections with 
the quantization conditions. Numerical TBA analysis shows a perfect match for the scal- 
ing functions obtained by the first method. Our results show that these two methods 
are complementary to each other. While the reflection amplitudes are confirmed by the 
numerical TBA analysis, the analytic structures of the TBA equations become clear only 
when the reflection amplitudes are introduced. 
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1 Introduction 



In the study of the two-dimensional quantum field theories near critical points, perturbed 
conformal field theory (CFT) approach has been quite successful HJ. Certain perturba- 
tions maintain the integrability structures for the models so that one can use the exact 
S'-matrices and the thermodynamic Bethe ansatz (TBA) methods to compute various 
physical quantities, in particular, the scaling function of the effective central charge as a 
function of the size of a system R 0. In the vicinity of the ultraviolet (UV) fixed point, 
the scaling functions of various models behave in different ways. While most common 
behaviour is the power law corrections of R, slow flows due to the l/(a — lni?) 2 correc- 
tions have been found for various affine Toda field theories with an unknown constant 
a. We show in this paper that to determine the UV behaviour completely one needs an 
independent method to compute the scaling function. 

The independent method we need has been first constructed in a remarkable paper 
by A. Zamolodchikov and Al. Zamolodchikov || where they introduce the "reflection 
amplitude" of the Liouville field theory (LFT) in terms of the correlation functions of 
the exponential operators and their dual fields. With the reflection amplitudes explicitly 
constructed from the structure constants of the LFT, they derived the scaling function of 
the ShG model as a function of a momentum which is conjugate to the bosonic zero-mode. 
Considered as an integrable perturbation of the LFT, the ShG model provides a confining 
potential well for the zero-mode so that the conjugate momentum should satisfy certain 
"quantization condition" which relates the momentum with R depending on the details 
of the reflection amplitudes. Numerical analysis shows a perfect agreement between the 
two scaling functions, confirming the new method. 

Our objectives in this paper is to compare the two methods both analytically and 
numerically in a deep UV region for the Bullough-Dodd (BD) model, another integrable 
perturbation of the LFT and to extend whole formalism to the super symmetric case. 
Formalism to study TBA equations analytically in the UV region has been first presented 
in |4[] for the ShG model by changing the nonlinear integral equation into infinite order 
differential equation and has been extended to affine Toda theories in ]§, 0. However, 
these works have considered only leading corrections of order l/(a — lni?) 2 . According to 
our analysis, the real interesting feature arises only when one takes into account higher 
order corrections where the quantization conditions found in the reflection amplitudes 
approach appear in the analysis of the TBA equations as a hidden structure. 

Comparison of these two methods can be used as a tool to check nonperturbative 
relations between masses of on-shell particles and dimensional parameters appearing in 
the actions. This is because the reflection amplitudes defined as off-shell quantities depend 
on the dimensional parameters while the TBA concerns only the on-shell quantities like the 
particle masses. The relation for the ShG model in |7| and generalizations to (fractionally) 
super symmetric theories in || can be tested against the TBA equations. 

Another interesting point happens when we analyze the supersymmetric sinh-Gordon 
(SShG) model using the two methods. As a perturbed super-LFT, there are two re- 
flection amplitudes corresponding to the Neveu-Schwarz (NS) and Ramond (R) sectors, 
respectively. Being different, they generate different scaling behaviours. Interestingly, the 
scaling function for the (R) sector has at most the power law corrections only. While 
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we have well-defined TBA equation for the (NS) sector, that for the (R) sector is not 
estabilished. We suggest the TBA of the (R) sector and provide the justifications based 
on the behaviour of the reflection amplitude. 

This paper is organized as follows. We introduce in sect. 2 the reflection amplitudes 
for the LFT and super-LFT and show that one can interpret the amplitudes as the 
quantum mechanical reflection amplitudes of the wave function. As an independent check, 
we solve the zero-mode quantum mechanical problem for the super-LFT to derive the 
reflection amplitudes showing that they are consistent with the quoted results. In sect. 3 
we consider the off-critical integrable models of the ShG, the BD and the SShG models 
and compute the scaling functions using the reflection amplitudes methods along with 
the quantization conditions. We perform the analytic computations in sect. 4 for the 
various TBA equations generalizing the leading order computation upto several higher 
orders enough for us to conclude that we can find nonperturbative equations identical 
to the quantization conditions. Numerical confirmations of our results are presented in 
sect. 5 and some relevent discussions are made in sect. 6. 

2 Reflection Amplitudes 

In this section, we introduce the reflection amplitudes for the LFT and super-LFT by 
quoting references. We interpret the reflection amplitudes as a reflection of quantum 
mechanical wave functional of zero-modes. 

2.1 Liouville Field Theory 

The LFT has been studied actively due to its relations to 2D quantum gravity and string 
theory and been shown that it enjoys all the properties as a CFT. The LFT action defined 
on a large disk T of radius R — > oo is 



where b is the dimensionless Liouville coupling constant and the scale parameter /x is 
usually called the cosmological constant and the background charge — Q at infinity is 



Al = — 




[ \{d a< pf + 47r/ie 2 ^j d 2 x + / M + 2Q 2 In R 



Q 



b + i/b. 



The LFT is a CFT with central charge 



1 + 6Q 2 



and the dimensions of exponential operators 



V a {x) = e 2a<t>(x) 



given by 



A a = a(Q - a). 
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Since Vq_ q and V a have the same dimension, Vq_ g is called "the reflection image" of V a 
and vice versa. These two operators are dual to each other. 

The n-point correlation function of the exponential fields defined as a functional inte- 
gral 

Q ai ,^a n (x 1 ,...,x n )= f V ai { Xl )...V an {x n )e' ALW D<j) 



can be determined as a formal expansion of the cosmological constant and computed 
exactly. The result shows an interesting relation between the correlation functions G ai ,...,a n 
and GQ~ai,...,a„- As an example, consider the three-point function which can be written 

as 

Ga 1 , a2 ,a 3 (x 1 ,x 2 ,x 3 ) = \xi2\ 2lA \x 23 \ 2ll \x2i\ 2l2 C(ai,a 2 ,a 3 ) 

where C(ai, a 2 , a 3 ) is the structure constant and 71 = A ai — A Q2 — A a3 and so on. There 
are several independent methods to compute the correlation functions; functional integral 



H [L0|, the canonical treatment pT| , and on-mass-shell conditional . 

The reflection a — > Q — a of each of the three operators introduces the Liouville 
reflection amplitude Sl(P) 

C(Q - oti, a 2 , ot 3 ) = C(ai, a 2 , a 3 )S L (iai - iQ/2) (1) 

where 

. . / „ 9 ^-2iP/bT(l + 2iP/b)T(l + 2iPb) 
SUP) = - ( W (f 2 )) ^ _ J^ {1 _ 2 J } (2) 

with 

7(x) = T(x)/r(l-x). 

By construction, the reflection amplitudes can be also defined from a two-point function 

(V a (z,z)V a (0,0)) = SL{%a ~^ ,2) with (V a (z,z)V Q - a (0,0)) = ^p-. 

Consider LFT on a cylinder of circumference 27r with the cartesian coordinates Xi, x 2 
where x 2 along the cylinder is defined as the imaginary time and X\ ~ x\ + 2tt is the space 
coordinate. The Hamiltonian acting in the space of states A of LFT 

^ = -^ + ^0 + ^0 

generates translations along the time x 2 . The space of states A is classified in the highest 
weight representations of Vir ® Vir 

A = ®pA P 

where a conformal class Ap contains a primary state vp = V a with 

Q 

a = ^ + iP 



and satisfies 



L„t;p = L n v P = at n > 

Lqv p = L v P = (Q 2 /A + P 2 )v P . 
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The primary state vp corresponding to the exponential operator V a becomes the lowest 
energy state and its descendants are generated by the action of L n and L n with n < on 
vp. Also the terminology of the "reflection image" becomes clear since the operator Vq- a 
is represented by the primary state V-p. Right and left generators L n and L n commute 
and therefore Ap has the structure of a direct product of right and left modules. 

In the LFT, one can reformulate the conformal structure in terms of the "zero-mode" 
of the Liouville field <f>(x) defined by 

o 2ir 

As 4>q — > — oo in the configuration space, one can neglect the exponential interaction term 
in the LFT action so that one can expand 4>(x) as a free massless field (z — x\ + 1x2) 



~? n I n n J 



(X) = (p 

where we defined the momentum conjugate to the zero- mode 0o as 

i d 



V 

and the oscillators satisfy 
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^0 



r r- - 1 - —x 

The Virasoro generators are given by 

L n = akdn-k + (27> + inQ)a n , n/0 

L = 2Ya- k a k + Q 2 /4 + P 2 , (3) 

fc>0 

and similarly for L„'s. The space of states is now represented as 

Aq = £ 2 (-oo < cj) < 00) <g> T (4) 

where £2 is the two-dimensional phase space spanned by 0o an d its conjugate momentum 
V and JF is the Fock space of the oscillators. 

Any state s E A can be represented by a wave functional in the 0o ~~ > — 00 

asymptotic limit. In particular, the wave functional for the primary state vp corresponds 
to 

^p[^i)] = (e 2iP0o + 5(P)e- 2iP ^)|O) as O -+ -00 (5) 

where S(P) is the reflection coefficient of the asymptotic wave functional. One can check 
that the wave functional of asymptotic form Eq.fJD has correct conformal dimension 
by acting L in Eq.(^). The coefficient S(P) should be the reflection amplitude Sl(P) 
introduced earlier since the wave functional ^ v _ p for the primary state v _p should be 
Sl(—P)^v p to be consistent with Eq.(|T]) along with 

S L (P)S L {-P) = 1. 
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In this framework, one can check the validity of the reflection amplitude by taking 
semiclassical limit b — > and using duality. Since P is of the order of 0(b), one can neglect 
the oscillators and keep only the zero-mode 0o so that the Hamiltonian is approximated 

as 

1 1 B 2 
#0 = -T7T-^ + 2^ e 2 ^. 
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The exact wave function of 0o for the Hamiltonian is well-known whose asymptotic form 
as 0o —* ~ o° is given by Eq.(f|) with 

,'7r/iV 2 ^ /6 r(l + 2iP/6) 



b 2 ) r (i - 2%p/b) ' 

It is straightforward to check this result is consistent with the non-perturbative reflection 
amplitude Eq.(@) perturbatively. 

2.2 super-Liouville Field Theory 

Now we extend above formalism to the N = 1 super-LFT whose lagrangian is given by 
With the background charge Q 

Q = b+l/b 

the central charge of the super-LFT is 

c 5L = ^(l + 2Q 2 ). 

The (NS) primary fields of the super-LFT are given by 

V a (z, z, 9, 9) = <f> a (z, z) + 9ip a (z, z) + 00 a (z, z) - 99(p a (z, z) 

with dimensions ^ 

A« = -a(Q - a), 

and the (R) fields by 

where is the 'twist field' with dimension 1/16 so that the dimension of the (R) fields 
are ^ 

&a = — + ~a(Q - a), 
lb 2 

The reflection amplitudes of the super-LFT defined from the structure constants have 
been derived from the structure constants in [13, 14].[] The reflection amplitudes for the 
(NS) fields are 

s ^ p) = -[t^{—)) r ( i-^)r(i-f)' (6) 



4 We quote the results after some minor corrections in such a way that they are consistent with both 
classical results and TBA. 
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and for the (R) fields 



s « <p) -[t<{—)) r (i-«)r(i (7) 
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The super-LFT is a super-CFT which satisfies the usual super- Virasoro algebra. The 
space of states for the super-LFT can be expressed by 

■Ao = £ 2 (-oo < 0o < 00,^0) ® F (8) 

where the fermionic zero-mode appears only for the (R) sector and T is the Fock space of 
bosonic and fermionic oscillators. The zero modes appear in the super- Virasora generator 
Lq and So of the (R) sector in such a way that L contains the square of the conjugate 
momentum V like Eq.([3]) and So acts non-trivially only on the twist field. 

The primary state vp can be also expressed by a wave functional ty Vp [(f>(xi)] whose 
asymptotic form is given similarly as Eq.(|5|). The amplitude S(P) is either Sns(P) or 
Sr(P) depending on the sector so that the wave functional *$? v _ p is given by S(— P)^ vp . 

One can also check the validity of this expression by taking the classical limit of h — ► 0. 
Since P is small of order of 0(b), one can neglect the oscillator part in Eq.(|8[) and study 
only the dynamics of zero-modes. In the (NS) sector, only bosonic zero-mode appears so 
that the Hamiltonian becomes 

which is essentially the same as that of the LFT, hence the reflection amplitude becomes 

Sm{p) ~ ~ [t) m-tp/by 

On the other hand, in the (R) sector, additional fermionic zero-mode is introduced in 
the hamiltonian by [|TJ] 

= - (^-X + Try&V^' + 2mfib 2 ^oe brpo . 
\d(f>oJ 

Since the fermionic zero-mode satisfies 

{■0o,^o} = O, ^0 = ^0 = 7; 



we can represent it by 



2 



1-1 i 



and the Hamiltonian becomes 



H + 
H 
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Among two eigen-spinors, the lower energy-state can be obtained as 



x K 1/2 _ iP/b {x) + K 1/2+iP/b {x)\ \ = Ho 
/ ' 



where K v [x) is the modified Bessel function. By taking the asymptotic limit 0o 
one can find the non-vanishing component is given by 

$ ~ e lP ^ + S R (P)e~ iP ^ 



-oo, 



with 



, (P) (w\-* p n\ + iPlh) 



These are consistent with the exact result Eq.fjTJ) in the b — > limit. 



3 Scaling Functions from Quantization Conditions 

We consider the ShG, the BD and the SShG models as integrable perturbations of the LFT 
and the super-LFT and compute the scaling functions by introducing the quantization 
conditions for the conjugate momentum in such a way that one can relate the momentum 
to the scaling parameter R through the reflection amplitudes. 



3.1 Perturbations of the LFT 

the ShG model 

We start by reviewing the analysis of P for the ShG model or A\ afline Toda field theory 
defined first on a circle of circumference R with periodic boundary condition. By rescaling 
the size to 2tt, one can write the action as 




where /x ~ [mass] 2+26 is the dimensional coupling constant with b the coupling constant. 

We are interested in the ground state energy E(R) or, more conveniently, the finite-size 
effective central charge 

a p 

c eS (R) = E(R) (10) 

7T 

in the ultraviolet limit R — *■ 0. Since we are interested in the ground-state energy, only the 
zero-mode contribution counts. So the corresponding effective central charge at R — > is 
determined mainly by P 

c eS (R) = 1 - 24P 2 + 0(R) (11) 

up to power corrections in R. 

For the ground state energy, one can consider only the zero-mode dynamics where 
the wave functional of <fio is confined in the potential barrier due to the ShG interaction 
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term. The ShG potential introduces a quantization condition for the momentum P which 
depends on the finite size R. As R — > 0, in particular, the wave functional is confined 
in the potential well where the potential vanishes in the most of the region and becomes 
nontrivial at 260o ~ ± ln/z(i?/27r) 2+2fc2 near the left and right edges. Near these edges of 
the potential well, the potential becomes that of the LFT and the wave functional will 
be reflected with the reflection amplitude of the LFT introduced earlier. Therefore, the 
quantization condition is given by 

(R/2iry 8iPQ Sl(P) = 1. 

In terms of the reflection phase 5l(P) defined by 

S L (P) = -e i5 ^ p \ 

the ground state momentum is qunatized as 

S L (P) =7r + 4PQln^. (12) 

Thus determined quantized momentum will give the scaling function c e g(R) in the UV 
region by Eq. (jTT|) . To see this explcitly, one can expand the reflection phase in the odd 
powers of P, 

5 L (P) = S^P + 5 3 (b)P 3 + 5 5 (b)P 5 + ... (13) 
where the coefficients can be obtained from the reflection amplitude Eq. (^) as follows: 

X fU\ 4 1 ^2 Ant l r (24^2 ) r V 1 + 2+2fe2 ) 

oAo) = - In b -4Qln — - A 1_ + 7 E 

4V7T 

1 fi 

s 3 (b) = yC(3)(fe 3 + r 3 ) 

64 

S 5 (b) = -C(5)(b 5 + b- 5 ) 
with Euler constant 7e- Now solving Eq.flTB"!) iteratively, we get 

c cS (R) = l+ C j + ^ + ^ + ... (14) 

where 

I = ^(6) - 4Q ln(i2/27r) (15) 
ci = -24tt 2 
C2 = -48n%(b) 
c 3 = -487r 5 (25 3 (6) 2 -vr5 5 (6)). 

The Gamma functions appear in 5\ due to the relation between the mass of the physical 
particle and the coupling constant /1 in the action 



71 fl 



m „ ( 1 



±yfH" V2 + 2&V r l 1 + 2 + 26 2 



2+26 2 
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the BD model 



(2) 

The BD model is an integrable field theory associated with A 2 affine Toda theory 
and can be regarded as an integrable perturbation of the LFT [ 12j . The action is given 
on a circle of circumference 2tc with periodic boundary condition; 



A 



BD 



where /i ~ [mass 



2+2b' 2 



2n 



2tx) 



2+2b 2 



e 2H + ^ { 



2+b 2 /2 



2tt J 



and // ~ [mass 

2\/3r 



m 



r(i + 



b 2 

6+3b 2 



6+36 2 



2+b2 / 2 are related to the mass of on-shell particle by 

2/x'vrr (l + ■ 



finT (1 + 6 2 ) 

r(-& 2 ) 



a+3b' 2 



r(-f) 



6+36^ 



This model possesses asymmetrical exponential potential terms compared with the 
ShG model. In the UV limit, the exponential potential becomes negligibly small except 
in the region where 0o goes to ±oo. This means that the BD model is again effectively 
described by the LFT. The scaling function of the central charge is given by the same 
Eqs.(TH]|) and (|TT|) in the ShG model. It is the quantization condition that makes the 
difference from the ShG model, due to the asymmetry of the potential well in the left and 
right edges. The conjugate momentum P is now quantized by the condition 

d . -4iP(Q+Q') 

" S L (P)S' L (P) = 1, 

— > 6/2 for Sl(P) given in Eq.(§) and 
Q' = 6/2 + 2/6. 



2ttJ 

where S' L (P) is obtained by substituting 6 

Q = 6+l/6, 
Using the phase shifts defined as 

S L {P) = -e i5 ^ p \ 



S' L {P) 



JS' L (P) 



the quantization condition becomes 

5(P) 

where 

m = \{h{P) 



n + AQP\n^- 

Ztc 



(16) 



Q'). 



The relation between P and R in Eq.(|T^) gives the scaling function continuous 
function of R, Eq . QUI) , with Q replaced by Q and <5's with <5's defined by power series 
expansion of the phase shift in P 



S(P) 



8 X P- 

6, 6 2 
-In — 

6 2 



5 3 P 3 + S 5 P 5 + ... 



2(Q + Q') 



In 



rnf (l + 



b 2 



6+3b 2 



•,3b 2 



2V3T 



+ 7e 



33 



32 



-C(5)(6 5 + - 



(17) 
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3.2 the SShG model 



Now we consider an integrable model obtained as a perturbation of the super-LFT, the 
SShG model. By rescaling the size to 2tt, one can express the action of the SShG model 
by 



A 



SShG 



■2-k 

dx2 I dxi 
R\ 1+b2 



07T Z7T 



+ 2ijjb' 2 ( tptp cosh(b(j)) + vr/iV 



R_ 

2^ 



2+2b 2 



[cosh(260) - 1] 



In the UV limit, the exponential potential becomes negligible except in the region where 
(f>Q goes to ±oo. This means that the SShG model is effectively described by the super- 
LFT as R — > 0. From the ground state energy for the primary state labelled by P, the 
effective central charge can be obtained by 

c cS (R) = ~12P 2 + 0(R) (NS) 
= -12P 2 + £>(#) (R). 

For the (NS) sector, P corresponding to the ground state is determined again by the 
quantization condition coming from the super-LFT reflection amplitudes: 

<WP) = 7r + 2QPln^, (18) 

Z7T 

where S^s (P) is the phase factor of (NS) reflection amplitudes. This quantization condi- 
tion can be solved iteratively by expanding <5ns(P) in powers of P, 



5f 



°3 



6™ 



NS D 3 



fNS D 5 



-2 In 



71 fX 



-7 



I CP) (» 3 + i 

-f«5>(* + * 



+ 7eQ 



(19) 



To decide the phases completely, one needs a relation between /i and m for the SShG 
model. This is given in || by 



7T 



2^7 



1 + 6 2 



yrfe 2 
1+6 2 



sin 



Trfc 2 
1+6 2 



1+fe 2 



(20) 



In terms of these coefficients, one can find the scaling function in a similar way as Eq. (|T 

We will consider the (R) sector in the next section since there is a fundamental differ 
ence from the (NS). 
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4 Perturbative TBA analysis in the UV region 



We compute the scaling functions analytically in the deep UV region by extending the 
methods developed in ||, [|, [J to higher orders. We find a hidden connection between the 
TBA and the quantization conditions arising from the reflection amplitudes. 

4.1 the ShG and the BD models 

The TBA equations for the ShG and the BD models are given by (r = mR) 

r cosh 9 = e(9) + / <p(0 - 9') In (l + e"^) (21) 

where the scaling function is expressed with the 'pseudo-energy' e(9) by 

CeS (R) = J J cosh 9 In (l + e~ e{0) ) d9 . (22) 

The model dependence comes from the kernel. The kernel of the ShG model is given by 

<p(9) = -<S> B (9), 

where 

_ ,„ s 4 sin 7ix cosh 9 . _ _ b 2 

$ x (6) = — — , with B 



cos 2tcx — cosh 29 ' 1 + b 2 

and for the BD model 

b 2 

¥>(0) = $2/3{9) + ®-B/3(0) + $(B-2)/a{8), with B 



l + b 2 /2' 



Following 0|, we express the effective central charge as 

6 f°° 

c e s(r) = — J r cosh 9L{9, r)dO. 

by defining 

L(9,r) = \n(l + e- t{e ' r) ). 
Fourier transform of the kernel <p(9) 



<p{k) 



/ d9ip(9)e ikd = 2tt \l + £ {-i) n $nk n \ 

J -°° { n=0 J 



rewrites Eq.(|2"ID as an infinite order differential equation 



rcoshtf + In (l - e- L{d ' r) ) = ^ <^ n L (n) 

n=0 



where 

L^ n \9,r) = (d/d9) n L(9,r). 
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We are going to concentrate on positive rapidity (9 > 0) since the TBA is even in 9. 
One can extend the solution to negative 9 using this symmetry. In the UV limit (r — > 0), 
since r cosh 9 can be approximated as re e /2, it is convenient to rescale 9 and define a new 
function 

L{9) = L{9 -ln(r/2),r) 
so that the TBA equation becomes 



+ In 1 



-L(fl) 



(23) 



n=0 



The central charge 



6 



Ccfrin 



7T^ Jln(r/2) 

becomes, after integrating by parts using -^e 9 = e' 

^ 2n-l 



e°L(9)d9 



71 



3 ' x 

~2 



]>> 2n £ (-l) fc+1 L( fc )(r')L( 2 - fc )(r') + ^ L(r')L(r') 



\n=l 



k=l 



6 



- — -C(l - e- i(r,) ) + -L(r') ln(l - e~ L ^) + e r ' L(r r 

7T Z \ 2 

where r' = In (r/2) and £ is the Rogers dilogarithm function. 

To solve the rescaled TBA, Eq. (|23D , we neglect the driving term e e and regard e~ L ^ 
as the same order of . (We are solving L{9) around the plateau region since c e fr is given 
in terms of L"(r')'s). The leading order of the TBA becomes the Liouville equation, 



whose solution is given by 



Lo(0)=ln 



fsm 2 {a(9-p)Y 



\ 2a 2 ip 2 J 
and the lowest order of the effective central charge is given by 



(24) 



CeSir) 



l + ^(^ 1 Hr') 2 -2e- L ^) + 



71- 



12 2 

1 o^2d + ■ ■ 

7X A 



a and (3 are integration constants, which are to be fixed by additional input. We first 
note that due to the symmetry property L(9,r) = L(—9,r) or 

L{9) =L(21n(r/2) -9), (25) 

an important relation appears between a and (3 

2a(P-ln(r/2)) = n7T, (26) 

where n is an arbitrary odd integer and is fixed as n — 1. Using this relation, we can 

"cos 2 (a(#-ln(r/2)))~ 



reexpress Eq. (|2~4l) as 



L (9) = In 



2a 2 (p 2 



(27) 
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As r — > 0, j3 remains finite while a — > 0. In addition, a should be small, a 2 (p 2 << 1 to 
make L > 0. Still, a and (3 are not completely fixed at this stage. We need the correct 
solution which vanishes as 9 — ► oo (tail part of L) to fix the integration constant. The 
lowest solution Lq{9) does not satisfy the correct condition, since we restrict the solution 
on the plateau: 9 may be restricted in the region 



7T 



< a{9 - In (r/2)) < - - y/2a 2 <p 2 

such that Lq(9) is positive and decreases as 9 increases. The restriction of the rapidity 
domain does not introduce much errors in the c e g\ the error is the order of 0(r), 

6r 



poo 

c eff ~ — / _ e 9 L(9)d9 ~ 0(r). 

7T J — yJKpo 



Since it is hard to get the complete analytic solution on the whole rapidity, we may resort 
to other physical solution to fix a and (3. Before doing this, we solve TBA on the plateau 
perturbatively. 

The rescaled TBA equation (^) can be solved by expanding in a series 

oo 

£ = E £«• 

ra=0 

Since Vq\Q) and e~ n ^ a ^ are of the same order, the expansion can be regarded as the 
expansion in a. We give the explicit solution of the TBA up to the order of a s . The 
differential equations for L 2 , L4, and L 6 are given as 

+ ^ = e~ L °L 2 - \e~^ 

+ + ^Lf = e~ L ° (L 4 - h(j + L 2 e- 2L ° - \e~^ 

^ 2 lf + $,Lf + (piLf + = e~ Lo (l 6 - L 2 L 4 + hi 

+ le- 2io ( L 4 - hi) + -L 2 e- 3Lo - h- 4L ° 



2 V 3 

These equations are iterative in the sense that one can find the solution for L 2n by 
inserting the solutions of previous differential equations, L 2 fc, k — 1, . . . n — 1. Since these 
equations are inhomogenious, the solutions are linear combinations of special solutions 
and solutions to the homogeneous equations 

U 



(p 2 L'i n = e~ L °L 2n . 

The solutions to the homogeneous equation are given by 

L 2n = Co(l + 2 tana;) + C\ tanx, x = a{9 — ln(r/2)). 

The second term c\ should vanish due to the symmetric property of L and the first term 
can be absorbed into L by redefining constant a. Therefore, it is enough to consider only 
special solutions. 
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The special solutions are surprisingly simple and given in terms of derviatives of L , 



a.(2,2)£o 2) + a 2a (2,o) 

: ( 4 ) , „2„. . f(2) 



L 4 = 0(4,4)^0 + a a (4,2)-^0 + a fl (4,0) 



^6 = a (6,6)^o 6) + ^ a (6,A)L'o > + a 4 a (6)2 )^ z; + a b a (6i0 ) 



("1) 



= (2) 



where a(«)'s are constants fixed by <p2nS. 



fl (2,2) 
0(2,0) 
0(4,4) 
°(4,2) 

a (4,0) 
«(6,6) 

«(6,4) 

fl (6,2) 

0(6,0) 



<^2 

4 

^2 



2(^4 



^2 



2592^| 



59$ - 900^2 + 4428$ - 2880<^ 2 

■x ,~2 i ino,»2 



,^2 



81^1 

5 ,2 



- 18^4^2 + 108(^4 - 90(^6^2 
2 



4 



y?4 



'<£>6 
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¥>2 



3<^2 



1 



259200^1 
1 



3240^1 



225(^1 
4 

45y3f 



281<^ - 7794^(^4 - 468504<^ - 69840<^ 6 
+712800^4^6 + 756^(131^ - 360<p 8 ) 
7fil - 198<^4 - 9288^1 

-1080^(^ 6 + 10800(^4^6 + 72^(31^2 _ 35 ^ g 
(pi - 24$£ 4 - 984^| - 190^(^ 6 

+1500(^2^4^6 + 16(^2(16(^4 - 35(^ g 



$ - 14c^ 4 - 264^ - 60^(^ 6 



+360^2^4^6 + 24<£ 2 (4y2 4 - 5(/9 



This solution shows some remarkable properties in connection with reflection ampli- 
tude. First, the relation between a and (3 in Eq. (^) is not changed after this higher order 
correction, since it is due to the reflection-symmetry property of the solution Eq. (|25|) and 
the higher order solutions are given in terms of derivatives of L . This symmetry rela- 
tion has the exact same analogue in the reflection amplitude, namely, the quantization 
condition Eq.([[|). Indeed, Eq. (|26D turns out to be the quantization condition once the 
integration constants a and (3 are identified with P and 5 by 



2PQ 
S(P) 



a 



2a/3 + 2a In tt = tt + 2a In (r/2vr) . 
By setting the mass m = 1, we will identify r with R. 
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Now one immediately notices that a has the role of P and (5 the phase of the reflec- 
tion amplitude and a and f3 have a new nonlinear relation in addition to the symmetry 
condition Eq. (|26|) , 



(3 — — In 7r 



2a 



ln7r H 

2a 



(i) 5i + (4 



5 3 + 



(2f 



Second, the higher order solution gives null contribution to c e s(r) up to this order, 



12 



c e &(r) = 1 ^y? 2 a + 0(a 



10 \ 



7T^ 



(29) 



= and explicit 
= -32a 6 . So, 
The same holds 



independent of the details of the kernel (p. We need only (p Q = (f2n+i 
values of L^'s at 9 = r': Lf \r') = -2a 2 , Lf{r') = -4a 4 , Lf\r 
Eq.(p9|) becomes exactly Eq.([l4]) since y5 2 = tc 2 /2Q 2 in the ShG model 
for the BD model by replacing Q with Q since (pi = n 2 /2Q 2 . 

The fact that higher order corrections vanish upto these orders makes it very plausible 
to conjecture that all the corrections indeed vanish so that c e g has only the lowest order 
contribution: 



c c s(r 



1 - 



12 



71^ 



; (f 2 a 



and consistent with the quantization condition arising in the reflection amplitudes. The 
scaling comes in through the relations between a and (3 only, Eqs. ( [26]) and (|28|). This is 
consistent with the reflection amplitude consideration. 



4.2 the SShG model 

The particle specturm of the SShG model is a doublet of mass degenerate which form 
a supermultiplet. Their S'-matrix has been obtained from the Yang-Baxter equation 
|17| which includes one unknown parameter. This parameter has been related to the 



coupling constant of the SShG model by interpreting the supermultiplet as bound states 
('breathers') of the solitons of the supersymmetric sine-Gordon model |L8] . 



Since the conventional TBA analysis gives the effective central charge in which only 
the lowest conformal dimension of the theory enter, the TBA equation will give only the 
(NS) result. Explicit derivation of the (NS) TBA based on the non-diagonal S'-matrix has 
been done in by diagonalizing the transfer matrix using the inversion relation. The 
scaling function can be expressed by 

Ceff (r) = J J cosh 6 In (l + e~ £l w ) d9 (30) 
where the pseudo-energies are the solution of the TBA equation, 

ei(0) = rcoshfl- / — y{6 - 6')hi[l + e'^'X 
J 2% 

e 2 (0) = -J^cp(e-9')Hl + e-^} (31) 
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where 7 = b 2 /(l + b 2 ) and the kernel is that of the ShG model. 

Analytic computations in the deep UV region can be done similarly. A little compli- 
cation arises due to the coupled TBA Eq. fl3ll). Defining 



L = ln 



1 + e 



-ei(0-ln (r/2),r) 



and M = In 



1 + e 



-e 2 (0-ln (r/2),r) 



we have 



-In 
-In 



1 - e^ (e) ] + M - L + ]T <^ 2 „M (2n) (#) 

71=1 
OO 

1 _ e -JfirWl + 1 _ m + ]T y 2n tt 2n \B) . 



n=l 



At plateau (6* — > —00) 
and at the edge {6 — ► 00) 



L(r') = M(r') , 

M(oo) = In 2. 



L(oo) =0, 
The effective central charge is given as 

Cjvs(r) = Co(r) + c fc (r) 

where 

ln(l - e~ L ) - ^ 2n L (2n) ] 
L - M - V2n{L {2n) - M (2n) )} L {1) d6 



coir J 



6 r°° 

Ti 2 Jr' 

6 



(32) 



Co(r) turns out to be the same as the central charge of the ShG model. This is because L 
vanishes at the edge and L(r') is described by the same TBA of the ShG model Eq. (p3|). 

Cfc(r) has the dominant contribution from the kink side of L(9) (non- vanishing part of 
L«). To express Ck in terms of L(r') we need two identities. One is given by integrating 
by part and substituting M with L using Eq. fl32|) 



ln(l - e 

T2n-1 



+ $^2n 



(-lY k+1) M (2n ' k \r')L (k \r') 



L k 



The other is given in terms of Rogers dilogarithmic function 



ln(l - e~™)M^de - -M(oo) 2 = c(- 

2 V 2 



Plugging these two identities into Eq. (|32|) we have 



00 g— nM(r') 



E 

n=l 



1 3 



7T" 



71, A.: 



rr 
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Combining with cq we have the scaling function of the SShG model, 



c(r) = 3_24^ + oWsjr) 



where we used the perturbative series results in the ShG model, a^s is the integration 
constant corresponding to a and (3ns to (3. Now recalling the reflection symmetry Eq. (|25|) 
and quantization condition of the reflection amplitude Eq. (|18D , 

2a NS {(3 NS -\n{r/2)) = n 
6 NS {p)-2QP]n^- = TT 

271 

we have the relation between the integration constants 

&NS = QP 

2a NS (3 NS = 5 NS (P) + 2QP\nir. 

With the help of (f>2 = tt 2 / (2Q 2 ) for the SShG model, the effective central charge is given 
by 

c NS = 3 --l2P 2 + 0(aZ). 

It is not clear how to derive the (R) sector TBA directly. Instead we can conjecture the 
TBA equation by considering slightly different version of the SShG model.0 By replacing 
the cosh super-potential with sinh one, we can still maintain the supersymmetry and 
integrability and can still consider the model as an integrable perturbation of the super- 
LFT. The particle spectrum, however, completely changes. Instead of mass degenerate 
of a boson and a fermion, we have only left- and right-moving massless fermionic modes 
where the boson becomes unstable and decays into two fermions. The S"-matrix between 
the two modes is the same as the ShG model. It is straightforward to write down the 
TBA equations since the scattering is diagonal. Now consider the (R) sector by imposing 



anti-periodic boundary condition. Following [16] for the Ising model, we can find the TBA 
of the (R) sector as 



ei (0) = rcosh#-^*ln(l- e - e2 ) (0) 
e 2 (0) = -<p * In (l - e~ ei ) (6) . 

The effective central charge becomes 

3r 



or f / 
Ceff(r) = — J cosh 6 In (l - e" ei 



iff) 



At the UV fixed point, one can express it in terms of Rogers dilogarithmic functions as 

e a e~ ea 



6 x . , 



In 1 - e~ ea - 



1 



6 

IT-" ' 



5 We thank Al. Zamolodchikov for suggesting this idea. 
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where variables 



are solutions of simple algebraic equations 



One can easily check 



X\ = 1 — X2, X2 = 1 — Xi 

V\ = 0, y 2 = 1. 



6 

c e ff = —z + £&) - C(x 2 )] = 



using a well-known identity 

£(x)+£(l-x) = £(l). 

In the UV region, there are only power corrections O(R). To see this, we rewrite the 
TBA in terms of H = In (1 — e~ e ) at the plateau (6 ~ r') by 

oo 

ln(l-e^) = # 2 + 5>2n^ 2n) 

n=l 

oo 

ln(l-e^) = H 1 + Y,V2nH{ 2n) . 

n=l 

For simplicity let us consider above equations with the assumption that H\ = H2 = H. 
By solving the leading order equation ln(l — e H ) = H (e H = 1/2), we can find higher 
order corrections around this. The next order correction satisfies the Hook's equation, 
(f>2^.H^ = —2 AH, whose general solution is given by AH = Asmuj(9 — (3) with u 2 = 
2/<f>2- Now this solution should satisfy the symmetry Eq. (pq), 

2u(/3 — In (r/2)) = (integer) x tt, 

It is impossible to satisfy this condition with finte u as r — > 0. Therefore, the correction 
AH should vanish and c e fj = + O(R). The same conclusion goes with the general case 

This result is equivalent to fixing the integer n appearing in the quantization condition 
to zero. The physical meaning becomes clear if one considers the P —>■ limit where 
Sr(P) — ► 1 comparing with Sns(P) — 1- While for the (NS) sector \l/p ~ 2iP(fi so 
that the quantum number n should be 1 as in Eq.(p!8|), the wave functional for the (R) 
sector becomes constant corresponding to n — 0. Therefore, the quantization condition 
becomes 

6 K (P) = 2Q P In 

Obvious solution is P = so that 



c cS (R) = + O(R). 
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In the b — > limit, one can verify this from the (R) sector zero-mode dynamics of the 
SShG model which is governed by the Hamiltonain 

H^ 1 = - ( Si- I + 47T 2 u 2 & 2 sinh 2 b(f) + Amub 2 ijj i> cosh&</> . 
\d(f>oJ 

This is a typical supersymmetric quantum mechanics problem and in general there ex- 
ists a zero-energy ground-state PD| if the supersymmetry is not broken. Explicitly, the 
wavefunction of the state is found to be 



^o(0o) 



g— 2n fi cosh b<j>o 





This state is normalizable and its energy is exactly zero. Thus at least in b — > limit, c e g 
is exactly zero regardless of r without any power correction. 



5 Numerical Analysis 

In this section, we solve TBA equations numerically and confirm the results in the previous 
sections. 

First we obtain the scaling function c e ff(r) of the BD model by solving numerically the 
TBA equation. Then, this result can be used to produce the reflection amplitude of LFT 
through the quantization condition which relates R(= r) and P in the following way. In 
r — > limit, one can neglect the 0(r) power correction in Eq. (pTTj) and define P(r) in the 
TBA framework as 

P = ^W- < 33 > 

The reflection phase S^ TBA \P) from TBA is now defined as the quantization condition, 
namely, 

5( TBA ) = vr + 4QPln— . (34) 
2n 

According to Eq. ([TT]) , 5^ TBA \P) should reproduce the Liouville phase 5(P) Eq.(|TjJ) up to 
exponentially small corrections in 1/P, 

5(P) =5 (TBA) + <D(e~^Q). 

In Table 1, we show the first three coefficients of <5( TBA ) in the expansion in powers of 
P obtained by numerical analysis of the BD model and compare with the corresponding 
^( LFT ) given by Eq. (|17D. We see an excellent agreement which confirms the validity of 
our approach. 

We have also plotted in Fig.l the scaling function c e ff(r) obtained from numerical anal- 
ysis of TBA equations and that from LFT reflection amplitudes (|TT|) , or equivalently, the 
full analytic evaluation of TBA equation (p9|) with 0(r) power corrections neglected, we 
find that they agree for r ~ 0.1 beyond which 0(r) power correction becomes important. 
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Figure 1: Plot of c cff for the BD and the SShG models. B = 0.5 for the BD model and 
B = 0.3 for the SShG model. 
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Table 1: First three coefficients of 5( TBA ) in the expansion in powers of P obtained by 
numerical analysis of the BD model in comparison with the corresponding <5( LFT ) given 
byEq.(|TJ. 

We can do the similar analysis for the (NS) sector of the SShG model. We define the 
momentum P in TBA framework as 

p =l/n (!-*<■■>)• 

and 5( TBA ) as 

5^ BA \P) = n + 2QP\n—. 

2n 

Table 2 shows the first three coefficients of £( TBA ) in the power expansion in P obtained 
by numerical analysis of the SShG model and the corresponding 5^ SLFT - ) given by Eq. (|T9|) 
supplemented with the // — m relation (p0[). Here again we see excellent agreements 
between the numerical and the analytical results. This result shows that the method 
based on the reflection amplitudes works perfectly in the presence of fermions like the 
SShG model. Especially the quantization condition for P holds perfectly even for b 
not small for which fermions interact nontrivially with bosons and the validity of the 
quantization condition may not be entirely clear at the first look. Thus our numerical 
result fully supports the analysis of TBA equations in sect. 4 and also the identification 
of undetermined constant j3 as that coming from SLFT reflection amplitude. Actually, 
from the plot of c e g(r) in Fig. 1, the approximation c e fj = 3/2 — 12P 2 is seen to work for 
larger region of r than that in the BD case. The supression of 0(r) correction might be 
related with the supersymmetric nature of the model. 



6 Conclusion 

In this paper we have analyzed the scaling functions of the ShG, the BD and the SShG 
models considered as integrable perturbations of the LFT and super-LFT using the re- 
flection amplitudes and conventional TBA analysis. Our main result is that the new 
method based on the reflection amplitudes is not only consistent with the TBA but also 
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Table 2: First three coefficients of <5( TBA ) in the expansion in powers of P obtained by 
numerical analysis of the SShG model in comparison with the corresponding 5( SLFT ) given 
byEq.flrjD. 
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necessary to completely understand the UV behaviour of the TBA. As our perturba- 
tive computations of the TBA show, there are no higher order corrections in the scaling 
functions except the leading correction which includes an unknown constant. To fix this 
constant, one should compare the scaling function obtained by the reflection amplitude, 
in particular the quantization conditions with the condition that the L function in the 
TBA analysis be symmetric. Indeed, we identified the hidden condition in the TBA which 
are equivalent to the quantization condition. Our analytic results are fully supported by 
numerical studies where we have correctly reconstructed the reflection amplitudes from 
the TBA. 

In addition, we find that the UV behaviours of the (NS) and (R) sectors in the SShG 
model are quite different. The scaling function of the (R) sector has only the power law 
perturbative corrections which vanish as b — ► 0. 

There are some interesting open problems which can be answered in a similar method 
used in this paper. First, our analysis of the ShG model, the A± affine Toda field theory, 
can be extended to those with higher rank by generalizing the one-dimensional quantum 
mechanical quantization condition to the higher- dimensional one JH]. 

Second, the new method to compute the scaling functions in terms of the reflection 
amplitudes can be extended to give the scaling behviour of the central charge of non- 
integrable models which can be expressed as perturbations of the LFT and super LFT. 
Previously, the scaling functions of the central charges can be computed only when the 
models are integrable so that one can use the TBA methods. 

Another problem is to generalize the quantization condition and corresponding quan- 
tum mechanical interpretation so that the scaling function obtained in this way makes 
sense for the whole range of the scale R. To do this, one may need non-perturbative 
expression for the power corrections. 

So far we have considered only the bulk integrable models. There has been much 
progress recently in the integrable models with boundaries. The main physical quantity 
which can be computed using boundary TBA is the boundary entropy which can also 
flow as the boundary scale varies. Using a similar logic presented in this paper, one may 
compute the scaling function of the boundary entropy using the reflection amplitudes of 
the boundary correlation functions for the boundary LFT |22[ . 

We hope our approach made in this paper can provide an alternative approach to 
compute the scaling functions in the two dimensional quantum field theories. 
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